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We derive the conservative part of the Lagrangian and the energy of a gravitation- 
ally bound two-body system at fourth post-Newtonian order, up to terms quadratic 
in the Newton constant. We also show that such terms are compatible with Lorentz 
invariance and we write an ansatz for the center-of mass position. The remaining 
terms carrying higher powers of the Newton constant are currently under investiga- 
tion. 
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I. INTRODUCTION 

The post-Newtonian approximation to the 2-body problem in General relativity repre- 
sents the common approach to study the bound object dynamics in the weak curvature, 
slow motion regime, see pQ and [2] for reviews. The interest in detailed analytical study of 
the gravitationally bound two-body problem has been revived by several concurrent factors. 

On the experimental side the era of gravitational wave astronomy is expected to start 
in a few years by the advent of the advanced gravitational wave detectors LIGO and Virgo 
[3J. Among other signals, these detectors will be sensitive to gravitational waves emitted 
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by coalescing binaries, making urgent to derive the equation of motion of binary black 
holes with the highest possible accuracy. The knowledge of accurate template waveform is 
needed for matching filtering techniques, commonly used in data analysis, whose output is 
particularly sensitive to the time varying phase of the signal which must be computed with 
0(1) precision [I]. 

On the theoretical side, since few years now stable numerical-relativity (NR) waveforms 
emitted in the last 0(10) orbits of a binary system (plus merger and ringdown to a final 
Kerr black hole) are available, see e.g. [5] for a recent review. As NR waveforms cannot be 
extended too early in the inspiral phase, during which binary constituents are far apart and 
non-relativistic, it appears natural to combine PN and NR results to construct complete 
waveform models: the highest possible accuracy on the analytical PN side is required in 
order to reduce the number of cycles of the NR waveforms needed to build complete, realiable 
waveforms [6]. 

On the theoretical analytical side, before the present paper, the Hamiltonian ruling the 
conservative dynamics of a gravitationally bound, spin-less binary system, has been com- 
puted at third PN order in [7] in the Arnowitt-Deser-Misner coordinates, calculation later 
confirmed by [81410"] in harmonic coordinates, and by [11], where a different resolution of the 
source singularity has been adopted. More recently, using the effective field theory (EFT) 
approach developed in [12], see [13] for a review, the Lagrangean of the gravitationally 
bound two body system has been re-derived at third PN order in [T4], with the use of an 
algorithm automatizing the decomposition of the problem into the sum of terms associated 
to Feynman diagrams. The application of the EFT methods to the PN approximation of 
General Relativity has also made possible the derivation of new results for spinning systems 
in both the conservative [15] and dissipative dynamics [16] . 

Following on our previous work [13], the main result of this paper is the computation 
of the two-body effective action for the conservative dynamics at 4PN order and up to 
terms quadratic in Newton's constant Gn, and thus represents a decisive step towards the 
computation of the full Hamiltonian dynamics at 4PN order. Using the virial relation 
v 2 ~ GnM/t, being r(v) the relative distance distance (velocity) of the binary constituents 
with M the total mass, the terms contributing to the 4PN order can be parametrized as 
G 5 N n v 2n with < n < 4. 

The remaining terms at 4PN not presented here involve 569 Feynman graphs whose 
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computation requires the evaluation of integrals among which the most complicated ones 
are few integrals analogous to 4-loop momentum integral in quantum field theory. The 
computation of the full Lagrangean at 4PN order is well under way and will be the subject 
of a future publication. 

We also note that the 4PN conservative contribution of the back-reaction force has been 
computed in [T7j and re-derived within EFT methods in [T5] . 

The paper is organized as follows. In sec. [IT] the fundamental gravity Lagrangean relevant 
for the two-body dynamics in the PN approximation is laid down and in sec. |III| we give a 
brief overview of the effective field theory methods for gravity, which are applied in the rest 
of the paper to the evaluation of the relevant Feyman graphs giving contribution to the 4PN 



order up G 2 N . In sec. IV standard techniques are applied to get rid of terms quadratic in the 
accelerations, in order to obtain a Lagrangean with terms at most linear in the acceleration. 



In sec. [V] the energy is derived and in sec. VI the Lorentz invariance is checked for and 
the center-of-mass position of the two-body system derived, always at 0(G 2 ) order. We 
summarize and conclude in sec. IVHI 



II. SHORT-SCALE LAGRANGIAN 

We perform our calculation in the EFT framework along the lines of [13]; in this section 
we resume our notations. The starting point is the action 

S = Seh + Sqf + S pp , (1) 

the first and third terms being, respectively, the usual Einstein-Hilbert action [25] and the 
worldline point particle action 

S EH = 2A 2 / d d+1 x^ R(g) , S pp = m t ! dr t = m, j yj-g^xfidx? d< ,(2) 

i=l,2 J i=l,2 ^ 

with 

A~ 2 = 32nG N L d ~ 3 . (3) 



As dimensional regularization will be needed in the computationm, L is a reference length 
scale which must cancel out from physical results and Gn is the standard 3+1-dimensional 
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Newton constant. As to the gauge fixing term Sgf, we follow [T9] : 



Sgf 



-A 



(4) 



with T M = r^p"^, which corresponds to the same harmonic gauge adopted in [TJ. Still 
following [IS], we adopt the standard Kaluza-Klein (KK) parametrization of the metric 
[20] (a somehow similar parametrization was first applied within the framework of a PN 
calculation in [2"T]): 



0/» = e 



-1 Aj-/A 
A . /A e -c d 0/A 7 .. _ AjAj/A 2 



(5) 



with 7y = 5jj + cry/A, Q = and i, j running over the (i spatial dimensions. In terms 



of the metric parametrization (|5]), each world-line coupling to the gravitational degrees of 



freedom <fi, Ai, <Tij reads 



S 



pp 



m / dr 



-m 



(6) 



and its Taylor expansion provides the various particle-gravity vertices of the EFT. 

Also the pure gravity sector Sbuik = Seh + Sgf can be explicitly written in terms of the 
KK variables; we report here only those terms which are needed for the full 4PN calculation 
(only a part of which is performed in the present work): 



Sbulk 3 / d d +W=7 {\ [(Va) 2 - 2(V^) 2 - (a 2 - 2(^) 2 ) e^ 



A 2 



(V</>) 2 - 2 e"^" 



+ 



2 



2 

e a + — 
A 

i 



FijA % Ai + AA(V-A) ) e^ - c d 0AV(/> 

2 A 2 



(7) 



+2c d ( 0VA - AV0] + -f ( A,rU + 2A k Tf j - 2A l T{ k ) - c d 



A 2 



The form ([7]) of the gravitational action [26] is in agreement with the one derived in [20] 



III. INTEGRATING OUT GRAVITY 



We have now all the ingredients to obtain a 2-body effective action S e ff with manifest 
power counting in Gjy and v at the desired (fourth) post-Newtonian order. This can be done 
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by integrating out the graviton fields from the full action derived above 

S e ff = J D(pDaijDA k exp[i(S E H + S GF + iS pp )) , (8) 

to obtain an effective action S e ff in terms of the particle positions and thier derivatives. 
As usual in field theory, the functional integration can be perturbatively expanded in terms 
of Feynman diagrams involving the gravitational degrees of freedom as internal lines |27j . 
regarded as dynamical fields emitted and absorbed by the point particles which are taken 
as non-dynamical sources. 

In order to allow manifest v scaling it is convenient to work with the space-Fourier 
transformed fields 

W£(t) = J d d x W a (t, x)e~ ip - x with W a = {0, A h . (9) 

The fields defined above are the fundamental variables in terms of which we are going to 
construct the Feynmann graphs; the action governing their dynamics can be found from 
eqns. ( 6|7 ). 



By looking in particular at the quadratic parts, one can explicitly write the propagators: 
P\W${ta)W$,(t b )] = ^P aa 5 ab (27r) d 5 d (p + p')V(p 2 ,t a ,t b )6(t a -t b ), (10) 
where P*+ = -± P A ^ = 5 ij: = - (fofy + S a 5 jk + (2 - c d )5 %] 5 kl ) and 

t a , t b ) = + ^ + (11) 

P 2 - dtA b P VP P J 

is the full relativistic propagator, which have been expanded as an instantaneaus nonrel- 

ativistic part plus insertion terms involving time derivatives (which after hitting the e tp ' x 

factors bring in extra v factors) and the infinite series is truncated at the some finite order 

in every Feynmann diagram. 

The next step is to to lay down and compute all the relevant Feynman diagrams; this 

has been done by means of a new version of the Mathematica code [22], including the use 

of the Feyncalc software [23], which has been used to reproduce the 3PN Lagrangian [Hj. 
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i : 
i : 
i \ 
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Figure 1: The three diagrams contributing at order Gn- The (j), A and a propagators are rep- 
resented respectively by blue dashed, red dotted and green solid lines. Propagators and vertices 
must be considered at appropriate v order to pick the 4PN contribution of these graphs. 

A. Order Gn 

The three relevant diagrams at this order are shown in fig. [TJ the typical graph contribu- 
tion A to the effective action ^ has the following form: 

A ~ G N m im2 L d - 3 j Yl i -^^e i ^-^V 1 V 2 , (12) 

where integration is performed over times tip an d the measure of momentum integration is 
d d p/(27r) d , n is an integer number with n < 4, x\p = Xipfap), and V\p = ^2^1,2(^1,2)] are 
the appropriate point particle vertices, expanded at the approriate order in i>i,2, which can 
be read from eq. ([6]). The integrations can be performed via standard Fourier formulae and 
do not lead to any divergency in dimensional regularization. 

We report below the result of the sum of the diagrams; symmetrization under exchange 
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of particles +(1 -H- 2) is understood throughout all this paper 
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The rationale for this separation is that such terms can be further reduced via the double 
zero trick (see e.g. [21]) which consists in writing, for instance 



(14) 



at - a 2 = {ax - eqi + eqi) ■ (a 2 - eq 2 + eq 2 ) 

= (ai - eqi) ■ (a 2 - eq 2 ) + eqi ■ a 2 + a x ■ eq 2 - eqi ■ eq 2 



(15) 



(ai - eqi) ■ [a 2 - eq 2 ) + (a x • a 2 J 



being egi 2 the value of a x 2 in terms of r, i> 1 and t> 2 as dictated by the equations of motion 



at the desired PN order. Since the first term on the last line of eq. 15 gives a vanishing 
contribution to the equations of motion, it can be dropped from the Lagrangian, thus leaving 
only terms linear in the accelerations, which however carry higher powers of Gat, as eq\ t2 = 
0(Gn). The standard form of the 3PN Lagrangian, as it is for instance reported in eq. (174) 
of p], is obtained via reiterated use of this trick. 



B. Order G 2 N 



The 23 diagrams contributing at 0{G 2 N ) are shown in figure [2j These diagrams are either 
factorizable into two single-graviton exchanges, or give a contribution A to the effective 



10 



/ / 

J. !_ 



/ '■ v '. \ \ ; \ \ 
/ . <•/ '/ '/ <• <• '/ . ■ ;• 
j. > •" / - ^ ■- ■• ^ / s ' .■ - _ ^.L ■ • 




Figure 2: The diagrams contributing at order G 2 N arranged in 4 lines. Diagrams in the n th line 
enter the dynamics at nPN order. Propagators and vertices must be considered at appropriate v 
order to pick the 4PN contribution of these diagrams. 

action of the type: 

A ~ G 2 N m\m 2 L 2 ^ £ f (td t J 2n « (zd t J 2n >> (td t2 ) 2n W la V lb V 2 

n a ,n b ,n 2 •'Mloi*16,*2 

/gi[(p— q)-a!i(tio)+q-a!i(tii)-p-a;2(*2)] 
q (p-q)2K + i) q 2K + i) p 2(n 2+ i) ^^ - ~ - h) , (16) 

with the rii < 3 and V^Ai = [p, q, dj la , dt lb , dt 2 ) being the bulk vertices describing the 
appropriate three-graviton interaction. 

These integrals can be reconduced to standard ones, eventually after some integration by 
parts, leading to the following total result: 
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.r2 



724 
~15 



133 

To" 1 



r2 



v[v r 2 3 + 



266 

\ 
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2T 



r 2 
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where terms involving more than one power of accelerations and their derivatives are: 
G 2 N m\m 2 
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268, _ 4211\ 

log r 

15 5 1800 



a\.a 2 
32 



28 , 377\ r 

y logr + — I ra x + 
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.r-2 



107 



30 



V>2 + 



g 2 2 1 12 g 2 



i T V 

+ 
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3 
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v 2 .a 1 v i .a 2 }■ + 



(d-3) ' 



We write separately the divergent terms: 
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+ a 2 .ai 
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with u = v\ — v 2 and f = \J 47re 7 £ (7 ~ .577216 is the Euler-Mascheroni constant). £^ 2 
is not relevant to determine the 4PN Lagrangian up to G 2 N because it gives rise to 0(G 3 N ) 
terms via the double zero trick. As to its divergent part C pole , it will have to be canceled by 
an appropriate worldline redefinition, analogously to what happens at 3PN [8], thus leaving 
a finite result in the d — > 3 limit. It is however not possible to determine the exact form 
of such redefinition at this point, as it must be further constrained by the knowledge of the 
poles occouring at higher powers of Gat. 
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IV. THE 4PN LAGRANGIAN LINEAR IN ACCELERATIONS UP TO G% 

Before we can write the Lagrangian linear in accelerations up to term of order G 2 Nl that 
is the 4PN extension of the standard formula reported in [T], we have to take into account 
of the manipulations that led to such standard formula. Indeed, to make it linear in the 
accelerations, the double zero trick have repeatedly used, and we must compute the effects 
of such transformation at 4PN. 



A. 2PN double zeros evaluated on 2PN equations of motion 



The following 2PN term quadratic in the accelerations 

PN G N mim 2 r 



C" 



16 



(15di ■ Q-2 — CLiCt^) 



(20) 



produces the 2PN and 3PN contributions correctly included in the standard 3PN Lagrangian, 
and the following 4PN term: 
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(21) 
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B. 3PN double zeros evaluated on 1PN equations of motion 



The following 3PN term quadratic in accelerations 
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96 
29 



23 



15 
32"' 



1 



32 
7 



7 



a2-cti - ^ a 2 a i («i + v 2 ) r + — rw.aaa! 



_ 5 55 o 29 r o 

— ra, V\.Vo-\ v-, H v\ H 1/, ^, 

32 1 16 32 1 32 1 16 1 



v r 

v\ v! 



+ ^ |( .'': 

47 
16 



9 

'i 

32 

47 



1 



v[v r 2 + 



r2 



Ojf .02 



41 

32 



13 

— \ 
32 



17 



-f1.a1W2.a2 



-■U2.a1f1.a2 



32 32 
generate (as well as several 3PN terms) the following contribution: 
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C. 4PN double zeroes coming from Gl graphs 



Finally, acting with the double zero trick on the Cq reported in eq. ( 14 ) gives 
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D. Summing up 

Summing up all the contributions, we have: 

r 4PN 7 ,,10 | riPN , MPN I /-o^S \ /qk\ 

= o^6 V i + C G N + C g% + °( G n) , (25) 
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where Cq^ = CffJ* — Cq n can be read directly in eq.(13), while C^ N takes into account 
all the contributions previously discussed, and reads: 
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The above formula, along with the term explicitly written in equation (13), is the main 
result of this paper. 
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V. ENERGY 



The following expression for the Energy is derived from the Lagrangian by means of a 
generalized Legendre transform : 

rpAPN 63 10 . G N m x m 2 J 525 8 1 4 1291 6 731 6 2 3 2 3 

= 256^ + 1 128 Vl + 32 ^ ~ 128" W + 128^ 2 + 21^ ^ 

373 4 , . 2 765 4 4 1511 4 2,02/ \2 2 , 3 / .3 2 
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6 /181 r 295 r \ r 4 /35 r2 35 r r 21 r2 
- ^ ( 128* + 128*) * + ^ U* + 32** + T* 
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32 i 2 32 1 2 32 1 2 320 2 y J v w; 
In the above expression, which is conserved up to order G 2 if one neglects gravitational 
radiation, accelerations appearing at 4PN have been systematically substituted by their 
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Newtonian value, and we included the post-Newtonian corrections of the analogous substi- 
tutions done in [TU] to write the 3PN Energy. 

VI. LORENTZ INVARIANCE AND CENTER-OF-MASS POSITION 

A. Boosts 

As discussed in [TOj, the post-Newtonian dynamics must be reminiscent of the symmetries 
of the fundamental theory. In practice, it must reflect the invariance under boosts, which in 
this framework are implemented by the following transformation 

<5£i,2 = -Vt + (V ■ x lt2 ) vi,2 + 0{V 2 ) , (28) 

where V 1 is the boost velocity. 

One must then require that the equations of motion are invariant under boosts, that is 
the Lagrangian variation must consist only of a time derivative plus possibly double zero 
terms: 

6C 4PN = y . ^pn + double zeroeg _ (29) 

As the boost transformation does not depend on Gjv (nor on d), one can to some extent 
check Lorentz invariance order by order in Newton's constant; we find indeed that the above 
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equation is compatible with Z APN = ^vfxi + Z^f + Z^ N + 0(G%), with 
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and 
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(31) 



(we remind that sum under particle exchange 1 -H- 2 is understood). 



More exactly, we find that eq.(29) holds up to order Gn if we assume that the following 
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are part of the double zero structure indicated at the r.h.s. of the same equation. Although 
we cannot be sure that this is indeed the case without knowing also the G 3 N part of the 
Lagrangian, it is a nontrivial fact that the only terms quadratic in the accelerations appear 
above; moreover, our assumption implies that such double zeroes must be accounted for 
when checking for the invariance at order G 2 N and this is exactly what we found. 

In its turn, Lorentz invariance at such order holds if we postulate the following 0{G 2 N ) 
double zero structure 
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whose actual existence will be checked only after determination of the 4PN G% dynamics 
(again, it is encouraging that only terms quadratic in accelerations are left). It is worth 
noticing that in deriving these results we had to take into account the 1PN corrections to 
the double zero terms which are already present at 3PN. 



21 



B. Center of mass position 



Still following [TU] , we can use the result of the previous subsection, and the fact of having 
a Lagrangian linear in acceleration, to determine the center-of-mass position at order G 2 N . 
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As for the Energy, accelerations appearing at 4PN have been systematically substituted by 
their Newtonian value, and we included the post-Newtonian corrections of the analogous 
substitutions done in [10] to derive 3PN the center-of-mass position. 



VII. CONCLUSIONS 



We have computed the the conservative dynamics of a binary system within the framework 
of the PN approximation to general relativity at fourth post-Newtonian order up to terms 
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quadratic in the Newton constant. By a sistematic use of the effective field theory methods 
for non-relativistic general relativity proposed by Goldberger and Rothstein it has been 
possible to automatize the computation and compute the effective Lagrangean by summing 
the contributions of several Feynman diagrams. We have also computed the energy function 
and verified for Lorentz invariance, giving a non trivial consistency check of our calculations. 
This work is the first ingredient for the complete determination of the 4PN dynamics. 
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T l - k is the connection of the purely spatial metric , Fij = Aj^ — Aij and indices must be 
raised and contracted via the (i-dimensional metric tensor 7, although in the last three lines of 
eq. [7] one can use just 5y as neglected terms are not needed at 4PN order; on the other hand 
all the spatial derivatives are meant to be simple (not covariant) ones and, when ambiguities 
might raise, gradients are always meant to act on contravariant fields (so that, for instance, 
V-l= jtiAij and Fl 3 = ^Vh)- 

As we focus on the conservative part of the dynamics, internal lines in Feynman diagrams 
correspond to the so-called potential gravitons and no gravitational raqdiation is emitted. 



